Paired states in spin-imbalanced atomic Fermi gases in one dimension 
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A growing expertise to engineer, manipulate and probe different cold-atom analogs of electronic 
condensed matter systems allows to probe properties of exotic pairing. We study paired states 
of spin-imbalanced ultracold atomic system of fermions with attractive short-range interactions 
in one-dimensional traps. Calculations are done using the Bethe Ansatz technique and the trap is 
incorporated into the solution via a local density approximation. The thermodynamic-Bethe- Ansatz 
equations are solved numerically and different local density profiles are calculated for zero and finite 
temperatures. A procedure to identify the homogeneous-system phase diagram using local density 
profiles in the trap is also proposed. Such scheme would be immediately useful for the experiments. 



Recently, ultracold-atom realizations of condensed- 
matter-physics models have attracted a lot of interest. 
A growing expertise in engineering, tuning and prob- 
ing atomic systems allows to study exceedingly compli- 
cated systems. One of the current challenges of con- 
densed matter physics is to understand the different ex- 
otic paired states that are realized when paired parti- 
cles have different chemical potentials. This disparity 
could be due to different physical effects such as differ- 
ent masses of the pairing particles or external magnetic 
filed and is thought to arise in systems such as quasi 
one-dimensional and two-dimensional organic supercon- 
ductors [l| and heavy-fermion materials [2] to name a 
few. It is argued that the zero-center-of-mass-momentum 
Cooper pairs are destabilized due the population im- 
balance and different, unconventional, paired states are 
proposed to occur 0, 0]. Among these exotic states, 
the FFLO (Fulde-Ferrell-Larkin-Ovchinnikov) state, in 
which pairs have nonzero center-of-mass momentum and 
superconductivity (superfluidity) coexists with nonzero 
polarization, is one of the candidates Recently, sev- 

eral experiments have studied the interplay of pairing and 
polarization in tree-dimensional (3d) spin-imbalanced su- 
perfluid clouds 0, Hj]. Ongoing experiments are explor- 
ing different pairing mechanisms in spin-imbalanced ul- 
tracold atomic systems of fermions in reduced dimen- 
sions [9], where the FFLO state is believed to have a 
large parameter regime of stability. The ultracold cloud 
of atoms is subjected to a 2d optical lattice, which de- 
fines an array of quasi-ld tubes, which can be regarded 
as independent if the intensity of the laser beams is large 
enough to suppress tunneling among them. Thus, to ex- 
plore pairing properties of Id systems of spin-imbalanced 
fermions is very relevant. There have been many differ- 
ent recent stu dies of this system. Bosonization 1(1 El]]. 



Bethe AnsatzEMl] , mean field [3], QMC OS El and 
DMRG lil-2lj methods have been used to study this 
problem and the agreement is that the partially polar- 
ized superfluid region (or "phase") is the analog of the 
FFLO state in higher dimensions (see Fig. [H explained 
below) . Most of the past analysis was done for zero tem- 
perature, except the QMC and mean field calculations. 



In addition, these approaches find the regime of strong 
interactions between particles very challenging, while ex- 
periments are done preferably in that regime. Therefore, 
there is a need for finite-temperature calculations, which 
do not suffer of the above limitations. 
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Figure 1: Zero-temperature phase digram of the uniform Id 
system, showing BCS-like (Id SC), FFLO-like (Id FFLO) and 
(fully polarized) normal (p = 100% curve) "phases". The areas 
marked "unstable" correspond to parameter regimes, which do 
not support equilibrium states. Local density trap profiles for 
total polarizations P = 4% (circles) and P = 48% (triangles) 
are superimposed on the diagram. 

In this paper, we address the problem using the 
Thermodynamic-Bethe- Ansatz (TBA) equations [221 [23]. 
which provide a full finite-T description and are well 
posed for arbitrary interaction strengths. We also pro- 
pose a scheme to reconstruct the "phase diagram" of a 
homogeneous Id attractive fermion cloud from experi- 
mental profiles of local particle density and polarization 
in a trap. We start from the Hamiltonian of a single Id 
tube in a harmonic trap 
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where m is the mass of the atoms, lo z is the trap frequency 

and 5lD = 1 - C al/ a± C = |C(l/2)|/>/2 ~ 
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1.0326] is the effective interaction strength between par- 
ticles, which can be tuned using a Feshbach resonance 
and made attractive (giD < 0) [2J|. a 3 D is the 3d scat- 
tering length and a± is the transverse oscillator length. 
In what follows, we shall measure lengths in units of the 
harmonic oscillator length (h/muiz) 1 ^ 2 and energies in 
units of hu z /2. The dimensionless Hamiltonian reads 



4A£)<J(*-2 i )+£\ 
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where A = — > o (attractive interactions). 

We will use the solution of the uniform system and take 
trap effects into account via the local density (Thomas- 
Fermi) approximation. This approximation is widely 
used in studying trapped gases, but its violation has been 
observed in highly elongated 3d clouds @| due to interface 
interaction effects 
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26j. In Id, one rather finds smooth 
crossovers and, for attractive fermions, comparison with 
QMC calculations [17J indicate that the local density ap- 
proximation gives good results for density profiles. Most 
features of local trap profiles are faithfully reproduced, 
but for the short wavelength Friedel oscillations, which 
are washed out at finite temperatures. In this framework, 

can be 
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the Bethe Ansatz solution of the model 
used locally. Eigenvalues and eigenstates of the system 
(of size L and with periodic boundary conditions) are de- 
termined from the roots of the Bethe- Ansatz equations: 
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where j = l,...,N and n = 1,...,M (M < N/2), with 
N the total number of particles and M the number 
of minority ones. {k{\ are particle quasi-momenta and 
solely determine energy eigenvalues, E — ^ = 

{A„} are (effective) spin rapidities and determine mag- 
netic properties of the system. Analyzing the roots of 
the Bethe Ansatz equations one can find that there are 
three classes of solutions ji^]: real fcj-s, representing un- 
paired particles; bound states of two particles with op- 
posite spins, represented by complex kj = Aj ± iA, with 
real spin rapidity Aj; and complex spin rapidities form- 
ing strings A (r)j = A (r) + iA(r + 1 - 2j), with j = 1, r 
and with real A^ r ' the center of the r-string. 

In the thermodynamic limit, L — > oo (keeping N/L 
and M/L constant), it is possible to locally define 
root and hole density distributions for unpaired parti- 
cles, bound states and n-strings: pu' l \k,z), p r b h \k,z) 

and psi' 1 ' (A, z) . The partition function is defined as 
Z = Tre- H ' T = e- F ' T , where £ = T(z) = e{z) - 
Ts{z) — n{z)n t {z) — hn s {z) is the free energy density, 
fi(z) = po — z 2 is the local chemical potential, h is the 



applied magnetic field, and s(z), n t (z) and n s (z) are 
local entropy, total-particle and spin densities, respec- 
tively. Then the equilibrium root and hole density dis- 
tributions are locally determined by minimizing the free 
energy while taking into account the Bethe- Ansatz equa- 
tions. The resulting TBA equations 13; 13| are 



fu - fu = G * /ft - G * f s i, 

fb - fb = 2[k 2 - A 2 - p{z)]/T + K x * f u + K 2 * f b , 

fsn — fsn — $n,lG * f u + G * (fsn+l + <W/sn-l)> 

lim (K n +1 * fsn ~ K n * fsn+l) = -2H/T, 

n—>-oo 

where f a = ln(l + p h J p r a ) and f a = ln(l + p r J 1 p h a ) are 
functions related to probability distributions, K n (k) = 

tt and '*' indicates the con- 
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volution of two functions (these notations are as in 
Ref. [13|). In terms of the /'s, T(z) reads 

T(z) = ~2^[ dkf u (k,z) - - J dkf b (k,z). (3) 

It is possible to straightforwardly determine densities 
of different thermodynamic variables from the free energy 
using x(z) = — g^^-j 1 , where x — n t , n s , s and x — p,h,T 
are conjugate thermodynamic variables. One can write 
equations for the different density distributions of all ex- 
tensive thermodynamic variables in a compact form 
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Psm + Psrn ~ "m.l^ * Pu + U * [Psm+1 + "m + Psm-lli 



lim {K n+ i * Plm -K n * pl^+i) = -6x,njn, 
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where the corresponding density distributions are de- 
fined as pY = (-1) 4 *.- ±dx(Tf y ) and p x y h = 
(-l) 1 + <5 f>= 2-d s (Tf y ). Notice that the above system of 
equations for the particle density (x = n t ) coincides with 
the continuum limit of the Bethe- Ansatz equations. Us- 
ing these definitions, local densities of thermodynamic 
observables in the trap are given by 



x(z) = / dkpY{k,z) + 2 / dkp x b r {k,z) 



(5) 



The integral equations are solved numerically for differ- 
ent temperatures and interaction strengths (cf. Ref. [29]); 
different density profiles are plotted in Fig. [2j The ref- 
erence chemical potential, /io, and the magnetic field, h, 
are fixed by the self-consistency relations J dzn t (z) = N 
and J dzn s (z) = PN . Depending on P, trap profiles 
show different regimes that reveal properties of the zero- 
temperature "phase digram" of a uniform system. For 
small polarizations (P — 4% plots), an FFLO-like state 
in the trap center coexists with a BCS-like state at 
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Figure 2: Trap density profiles for N = 200, T = 0, 0.06, 0.12[/i|F] (dashed, solid and dotted, respectively) and different total 
polarization P. The first column shows spin-up (rif, solid), spin-down (nj., dashed) particle densities and their difference (n s , 
dotted) for A = 25 and T = only. The second column shows the local particle densities (nt) for A = 25, 100 and different 
temperatures. The third column shows the respective local polarizations (p). The n t and p plots for A = 25 are vertically 
offset for clarity. 



the edges; while for large polarizations (P = 26% and 
P = 48% plots), the FFLO-like state in the trap center- 
coexists with a fully polarized normal state at the edges. 
Close to the "critical point" (P = 15% plots), the entire 
trap is in the FFLO-like state. The interface between 
different "phases" in the trap is marked by a kink feature 
in the local particle density (or polarization) profiles at 
zero temperature, which is smeared out by thermal fluc- 
tuations. From Fig. [2] we see that while for cases with 
large total polarization the phase boundary is still visible 
at T = 0.06 \p^f] (here /ijf is the noninteracting Fermi 
energy for the majority species), for small polarizations 
temperature should be even smaller to observe the kink. 
From the local polarization profiles, we see that in the 
normal state temperature induces bound pairs, while in 
the fully paired state unpaired particles appear. 



We propose here a scheme to identify the low- 
temperature "phase diagram" of the uniform system from 
trap profiles of the local particle density and polarization. 
To achieve that aim we represent the phase diagram in 
variables that can be measured directly in the experi- 
ments. We argue that nt(z) and the majority-spin local 
chemical potential ^]{z) = fi(z) + h can be straightfor- 
wardly measured and, thus, a "phase diagram" in these 
variables should be directly useful to analyze experimen- 
tal findings. In these variables, the "phase diagram" of 
the uniform system (see Fig. [1} shows three different 
states: the BCS-like fully paired state with zero polar- 
ization (IdSC), the FFLO-like state with non-zero polar- 
ization (ldFFLO) and the fully polarized normal state 
which has collapsed into a curve (cf. Fig. 1 in Ref. [12]) • 
Plotting the local particle density and polarization pro- 
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files as functions of /u-f will reconstruct the uniform- 
system "phase diagram" and probe into the ldFFLO 
state. To elucidate the scheme we consider a trap profile 
with the total polarization P — 48% (triangles on Fig.[Tj. 
In this regime, the ldFFLO "phase" in the trap center 
coexists with the fully polarized region at the edges (see 
Fig. [21 . In the framework of the local density approxima- 
tion, the radius of the cloud is given by R = V^o + h and 
thus one can easily find ^\{z) + h = R 2 — z 2 . Thus, the 
trap profile can be uniquely placed on the "phase digram" 
(see Fig. [TJ. The above identification is valid for clouds 
with fully polarized wings and plotting many profiles of 
these type allows to probe into the ldFFLO. For small to- 
tal polarizations, the ldFFLO in the trap center coexists 
with the IdSC regions at the edges (see Fig. [2]). In this 
case, the radius of the cloud is given by R = \/A 2 + /i 
and thus one can easily find ^(z) + h = R 2 + h — A 2 — z 2 . 
It is again possible to place a trap profile on the "phase 
digram" (see trap profile for P = 4% in Fig. [TJ. Due to 
the nonzero temperature and measurement noise, it will 
be challenging to precisely identify the zero polarization 
curve (ldFFLO-ldSC boundary), but most features of 
the "phase diagram" (including inner ldFFLO isopolar- 
ization lines) are, nevertheless, within the reach of the 
upcoming experiments. 

To summarize, we have studied different pairing states 
in spin-imbalanced ultracold atomic clouds of fermions in 
Id. We have combined the Bethe Ansatz method and the 
local density approximation to calculate trap profiles of 
different observables. Our full finite-temperature calcula- 
tions shed light on the parameter regimes for total den- 
sity, polarization and temperature that are required to 
observe exotic states such as the ldFFLO. We also pro- 
posed a way to unambiguously reconstruct the "phase" 
diagram of the uniform system from experimental local 
density profiles of the trapped clouds. Our calculational 
scheme could also be used to compute asymptotic mo- 
mentum distributions that could be measured after si- 
multaneously releasing the trap and switching off the 
interactions in time-of- flight experiments. This would 
be complementary to the in-situ pair-momentum corre- 
lators discussed already in the literature [3-[i(l which 
are deemed to provide extra signatures of FFLO physics. 
A similar analysis can be also performed for repulsive 
interactions between particles, where the coexistence of 
different Luttinger-liquid regimes has been proposed [30l | . 
The above approach is well suited to explore interesting 
crossover effects, which are challenging to capture with 
other methods. Such a system would also be within the 
reach of the present generation of experiments. 
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